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ABSTRACT 

A statistical method is presented for determining the velocity field in the 
immediate vicinity of groups of galaxies using only positional and redshift in- 
formation with the goal of studying the perturbation of the Hubble flow around 
groups more distant than the Local Group. The velocities are assumed to obey 
a Hubble-like expansion law, i.e. V = H exp R where the expansion rate H exp is to 
be determined. The method is applied to a large, representative group catalog 
and evidence is found for a sub-Hubble expansion rate within two well defined 
radii beyond the virial radii of the groups. This result is consistent with that 
of Teerikorpi et al. (2008) who found a similar expansion law around 3 nearby 
groups and extends it to a more representative volume of space. 

Subject headings: galaxies:groups 

1. Introduction 

Teerikorpi et al. (2008) (hereafter T08) have recently shown that galaxies in the im- 
mediate vicinity of 3 nearby groups of galaxies exhibit a Hubble-like expansion law at a 
slightly sub-Hubble expansion rate. The authors had at their disposal measured velocities 
and distances in order to arrive at this result. T08 interpret their result in terms of the 
effects of dark energy. The arguments for the dark energy interpretation can be found in a 
series of papers starting with Chernin et al. (2000). Briefly, when a density perturbation 
which is to become a galaxy group collapses, it leaves behind an empty annulus between the 
original turn-around radius designated here Res, an d its post-collapse zero-velocity radius 
which lies between the virial radius and R v . R v is the radius where the gravity force of the 
group mass is balanced by the antigravity force of the vacuum energy. Should a galaxy find 
itself in the region between R v and Res either as a result of the chaos of virialization or 
by later ejection from the group, it may be accelerated by the effects of the vacuum energy. 
Since as these authors argue, the radial motion of the particles within this shell should be 
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nearly the same as Newtonian motion its expansion velocity can be described by the vacuum 
Hubble constant H v = y/U\H . Maccio et al.(2005) found from N-body simulations that 
the observed coldness of the local Hubble flow around the local group is consistent with the 
dark energy interpretation. 

However recently Hoffman et al. (2008) and Martinez- Vaquero et al. (2009) using 
sophisticated cosmological simulations claim that dark energy should have no effect on the 
local dynamics. In particular, these authors showed that a similar local Hubble flow was also 
found in an open cosmological model making the need for a A model unnecessary locally. 
Further, Sandage (1986) has shown analytically that sub-Hubble flows can exist around 
groups in cosmologies without a A term though of a different form from that found by T08. 
It is clearly of interest to confirm and extend the T08 result. In order to do so a new method 
for finding the velocity field is required since with increasing distances the only observational 
data available are positions and redshifts both for the groups and the field galaxies. The 
goal here is to devise such a method and to apply it to an existing group catalog in order to 
determine if the perturbations to the Hubble expansion law found by T08 apply to groups 
in a larger and more representative volume of space. 



2. Fitting the Velocity Field Around Galaxy Groups with a Hubble-like 

Expansion Law 

Determining the expansion velocity in terms of Hubble-like law (V = H exp R) requires 
knowledge of the velocity with respect to the group center in the direction of expansion V, 
and the distance from the group center R. To demonstrate the essence of the method, assume 
for the moment that two projections of this expansion law are available (i.e. AV = Vcos{9) 
and AR = Rsin{9) in Fig f ). Let y' — V cos{6) and x' = H exp R sin{6). Then 

y'jx' = -^—cos(9)/sin(9) (I) 

Now take the average of both sides of this equation assuming that the number of groups is 
large and that the surrounding field galaxies are randomly distributed about these groups 
and denote the mean of the angular term as 

I n = / cos{9) / sin{9)du (2) 

J sphere 

after rearranging terms H exp is varied until the ratio 



(y'/x')/I n = 1 = V/H exp R (3) 
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and the true expansion rate is obtained. In practice, instead of V cos{6) and R sin{9) which 
are not known, we do have access to the components V cos((3) and R sin((3) (Fig.l) so an 
intermediate step is required and is described below. 

The observational data available consists of the angular positions and redshifts of the 
group centers and the field galaxies (z g and zf). For a given cosmology, co-moving distances 
d(z) can be calculated from d{z) = c/H Q dz'/((l -fi A )(l + z') 3 + tt A ) 1/2 . Note that £l k = 
is assumed allowing the use of Euclidean geometry. Further the redshifts of all the data to 
be analysed are < 0.1 so that the changes due to different assumed cosmologies are small. 
From Fig.l with d g now available 

r p — R sin(/3) = d g sin(a) (4) 

and 

dper = d g cos(a) (5) 

Numerically inverting the above expression for co-moving distance allows the calculation of 
z per , the redshift at d per . 

Now consider the line of sight projection of the expansion law AVd where 

AV D = V COS(P) ~ C (Zf - Z per ) (6) 

and define the ratio y/x as follows 

y/x = AV D /(r p H exp ) = (V / H exp R)cos{f3) / sin{[5) (7) 

Rotate the line-of-sight about the observer by + a and define a new ratio in terms of 9 = a+(3. 
Since V and R remain unchanged by this operation we obtain 

y'/x' =| (—sin(a) + (y / x)cos(a)) / (cos(a) + (y/x)sin(a)) \— (V/H exp R) \ cos (9) / ' sin{6) \ 

(8) 

Note that by taking the absolute value in (8) the integral over the sphere becomes by sym- 
metry 

j = 27rf* /2 \(cos(9)/stn(9))\szn(9)d9 = ^ 
2tt j*' 2 sm{9)d9 

Hence 

V/H exp R = VJx 1 (10) 

where the true value of H exp is that which makes the right hand side of (10) equal to unity. 
For convenience we shall refer to the statistic represented by (10) as Y(H exp ). Lastly from 
Fig 1 and equation (7) 

R(H exp ) = r p /sin(/3) = r p /sin(cor l (AV D /(r p H exp ))) (11) 
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where H exp is the value obtained above. 

The above analysis applies to any one region undergoing Hubble-like expansion at a rate 
H eX p- 

The region with the sub-Hubble expansion rate (H v ) found by T08 lies between R v and 
Res hi Fig.l. (R v , Res an d H v are defined in the next section). Beyond Res the expansion 
is assumed to return to the global Hubble rate H Q (H(z)). (At higher z the global rate 
becomes H(z) = H ((l — ^a)(1 + z) 3 + fi A )^ 2 )- Similarly, by (11) in each of these regions 
we then have R(H V ) and R(H(z g )). 

With two expanding regions, the four key quantities in the analysis are the mean of 
the ratio y'/x' (equ'n 10) evaluated with H v and H(z g ) replacing H exp designated here 
Y(H V ) and Y(H(z g )) and the normalized distances R(H V )/R V and R(H(z g ))/R v . These four 
quantities are calculated for every field galaxy. In order to confirm the T08 result, Y(H V ) 
should be unity when 1 < R(H V )/R V < Res/Rv for the particular value of H exp = H v and 
Y(H(z g )) should be unity when R(H(z g ))/R v becomes larger than R E s/R v - Both of the 
above conditions rely on the assumption that the galaxies are randomly distributed about 
the group centers for the sample as a whole. Note also that because we are dealing with 
ratios, the result is independent of the actual value assumed for H D . 

In order to validate the above procedure a simple numerical simulation was performed. 
In this simulation the same group properties i.e. the positions, redshifts and masses (from 
which Ry could be calculated) of groups used in the analysis of the observations were used as 
centers around which 'field galaxies' were randomly distributed (both in angle 9 and group- 
centric distance R). Around each of the 12620 group centers are 300 randomly chosen field 
galaxies. For a given cosmology, line-of-sight redshifts were assigned (either as a result of 
global Hubble expansion or sub- Hubble expansion about the group center) depending on 
whether R/R v was greater than or less than 2. The simulation is clearly idealized to the 
extent that each group and its retinue of randomly distributed galaxies was treated as an 
isolated system. The data were then analysed (as described above) using the same routines 
that were used for the observations. Here the values of R plotted on the abscissa come from 
the simulation. Shown in Fig. 2 are the results for an assumed (fi m , f^v, Qk) (0.3,0.7 ,0.0) 
cosmology. The red points represent Y(H V ) and R/R v . The black points were calculated 
assuming H exp = H(z g ) and plotted as Y(H(z g )) versus R/R v . Note the characteristic 
step in the distribution of both sets of points at the transition point (R/R v = 2) and the 
continuity of red and black points at V/H exp R = 1. 
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3. Analysis of the Observations 
3.1. The Data 

The Tago et al. (2006) group catalogue of galaxies from the 2DF redshift survey (Colless 
et al. 2001, 2003) consists of 3 tables-one enumerates the groups and their properties, the 
second lists the group members while the third gives the co-moving distances of the galaxies 
which do not belong to any group. This catalogue was constructed using a friends-of friends 
algorithm. As such any individual group is subject to contamination problems (e.g. Niemi 
& Valtonen 2009). In the application discussed here, the groups are used statistically so that 
any contamination effects will affect the determination of both sub-Hubble and the global 
Hubble rate in the same way. For this investigation we limited our selection of groups to 
those lying within 300/h Mpc. This sample contains 12620 groups. The virial theorem (i.e. 
M = (tt /2)r viria i3a 2 /G) where r viria i and a are given in the first table of the group catalog 
was used to determine the mass of each of the groups. These masses were then used to 
calculate values of R v from equation (12) below. 

For a given group each field galaxy subtends an angle a between it and the group center 
and each lies at the tabulated co-moving distance computed for a (fi m , £Ia, ^fc) (0.3,0.7,0.0) 
with H D = 100/i and h = 1 cosmology. These data are used to recover the individual 
redshifts (z g and z/) which the authors had corrected for the motion relative to the CMB. 
The redshifts are then used to recompute the co-moving distances d g and df (now assuming 
h=0.7) for the different trial cosmologies considered below. 



3.2. Some Definitions 

For the purpose of comparing the velocity field derived here with that of T08, we define 
the following analogous quantities. 

Let M represent the mass of an individual group and define the two radii R v , and Res 
as follows. 

M/(4tt^/3) = 2p v = 2Q APc (12) 

and 

M/(4ttR es /3) = Pm = (l- n A ) Pc (13) 

here p c is the critical density and is equal to 3H^/8nG. R v is the radius where the gravity 
force of mass M is balanced by the antigravity force of the dark energy 0. Res is the 



The factor of 2 in equation (12) arises because in General Relativity the effective density includes a 
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radius within which the the mass now in the group originated and beyond which the velocity 
field approaches the global Hubble flow (Teerikorpi & Chernin, 2010 hereafter T10). Note 
that Res/Rv = (2^a/(1 — f^)) 1 / 3 and is greater than one as long as f^A > 1/3. The 
region within Res is referred to by T08 as an Einstein-Straus vacuole. Also recall that 
M/(4irR 3 virial /3) ~ 200p c so that R v /R mna i ~ 5.1. Rearranging (12) 

R v = (GM/(H 2 o n A ))^ (14) 

Finally we parameterize the unknown sub-Hubble expansion rate as 

H v = a/^aX (15) 

where following T08 we have neglected higher order terms (Chernin et al., 2000). This 
expression (referred to as the vacuum Hubble constant) also follows from the first Friedmann 
equation with fl m and Qk set to zero and is independent of redshift. 

Before presenting the observations, we discuss some limitations of the method imposed 
by the observations themselves. Velocity components along the line of sight beyond the 
projection of the actual expansion velocity will bias our statistic to be larger than one. 
Two contributors are redshift errors and peculiar motions (a 'noisy' Hubble flow). Redshift 
errors are ~ 85 kmsec^ 1 (Colless et al. 2001). The imprecise location of the group center will 
generally bias our statistic to be less than one. If both the inner and global expansion laws are 
affected in the same way the severity of the above biases can be monitored by evaluating the 
statistic for galaxies undergoing global Hubble expansion. Another point to consider is that 
groups that are close together can have the same field galaxies in common. However, just as 
galaxies participating in the Hubble flow are doing so irrespective of the co-moving center, the 
same argument applies to galaxies assumed to be expanding within the vacuum dominated 
region. For this reason, we choose galaxies within each region 1 < R/R v < R E s/R v (vacuum 
dominated region) and R/R v > Res/Rv (global Hubble expansion region) and treat them 
as separate samples. 

In order to determine the dividing line between sub-Hubble and global expansion, solu- 
tions were obtained for the data broken up into small groups with similar values of R v . It was 
found that the distinct step pattern seen in the simulation (Fig. 2) while visible for the whole 
sample was most clearly defined when using only those groups with 2.6 < R v {h^Mpc) < 3.3 
for Q\ = 0.7. The corresponding range in group mass is ~ 1 — 2 x 10 13 h~ 1 M & . The results 
are shown in the top two panels of Fig. 3. Similar results were found for other trial values 



contribution from the pressure as well i.e. p e ff = p + 3P/ c 2 and P = 



-PvC 2 . 
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of Qa after allowing for the JIa dependence of R v and are shown in Fig. 5. The bottom 
panel of Fig. 3 shows the results of assuming an open cosmology model (OCDM) with 
Vt rn = 0.3, Q\ = 0, Qk — 0.7 for the same groups and with the above values of R v used here 
in order to allow a comparison of results between the two cosmological models. Technically 
R v = oo for the open model. Common to all 3 panels, the point-to-point scatter appears 
smaller below the break than it does above. The dispersion of the distribution of Y(H(z g )) 
for 117896 group-galaxy pairs above the break (middle panel) is twice that of the dispersion 
of Y(H V ) for the 9894 group-galaxy pairs below the break (top panel). Comparable numbers 
apply to the other models shown in Fig. 5. If the galaxies below the break originated from 
within the group they will have the group peculiar motion imprinted on them and they will 
exhibit a 'cooler' flow. Beyond the break the galaxies are uncoupled from the groups and 
show the full effects of peculiar motion. This applies here to both the OCDM and ACDM 
model because the break appears in the same place. The bottom 2 panels show that the 
global Hubble rate cannot account for expansion velocities of the galaxies in the region below 
the break. The ACDM model can provide the appropriate sub-Hubble rate (top panel). 

Fig. 3 shows that the dividing line between the two velocity regimes is reasonably 
sharp and occurs at R/R v = 2.0 ± 0.1. We identify this particular value of R as Res- 
(Note from equ'ns (12) and (13) the ratio Res/Rv is predicted to be 1.7 when f2 A = 0.7 
and 2 when f2 A = 0.8 (T10). A possible explanation for this small discrepancy comes from 
the observations of galaxy clusters that a substantial amount of mass (up to 100%) exists 
outside the virial radius but within the post- collapse zero-velocity radius which is below R v 
(e.g. Rines & Diaferio 2006). Allowing for this possibility would increase R v by 2 1 / 3 and 
the predicted ratio would be recovered). Another possible contributor to the discrepancy 
in Res/Rv is our neglect of the 'lost-gravity' term due to dark energy in the calculation 
of the virial mass (Chernin et al. 2009). It is instructive to divide the data at this point 
and to plot the statistic (equ'n 10) versus R v itself. The results are shown in Fig. 4 for 
the case fl\ = 0.7. The top panel shows results for R/R v < 2 for all 12577 groups with 
R v < 5.7hj^Mpc. The middle panel shows results for the same sample for R/R v > 2. 
As before the red points represent the statistic calculated with H exp — H v — 0.8AH o (i.e. 
f2 A = 0.7) while the black points were obtained with H exp = H(z g ). The dotted lines enclose 
the region used to construct Figs. 3 & 5. The region at R v < 2hj^Mpc never reaches 
unity while the data become very noisy at R v > 3.3hjQMpc and the break becomes less well 
defined. It is interesting that the N-body simulations of Wang et al. (2009) show (their Fig. 
6) that for host halo masses greater than 10 l3 h~ 1 M Q ejected sub-halos fall back while for host 
halo masses lower than this ejected sub-halos are likely to reach the region beyond R v . The 
data in this higher mass range are too sparse to determine whether this same phenomenon 
is responsible for the break becoming less well defined. Plots for higher and lower values 
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of Q\ are similar with the statistic among the red points between the dashed lines in the 
R/R v < 2 plots asymptoting at lower and higher values. The results for R/R v > 2 (black 
points) for a range in vacuum energy are very similar to those shown in the middle panel. It 
is not clear what is responsible for the fall-off below R v ~ 2h^Mpc in these plots especially 
given that T08 found their result at R v ~ 1.3/ifo 1 Mpc. A possible explanation comes from 
the observation that the galaxies that are supposed to be expanding at the Hubble rate 
(those plotted in the middle panel) also show large deviations from unity in the same low 
range of R v suggesting that the field galaxies surrounding these low mass groups may not 
be uniformly distributed. 

A sub-sample of more isolated groups was analysed in order to measure the effect of 
group clustering on the result. This sub-sample included only those 3894 groups of the 
original 12620 without a neighbor within Ah^Mpc. The results are shown in the bottom 
panel of Fig. 4 and can be compared to the top panel. The previous results shown in 
Figs 3 & 5 remain unchanged but with larger uncertainty. In addition the fall-off below 
R v ~ 2hjQMpc noted above is still present. This region should show the most change if the 
effect of clustering is a dominating influence. Increasing the isolation criterion above leaves 
too few groups and very poor statistics. 

As a further check on the results another sub-sample consisting of only those 1835 groups 
with more than 5 members was analysed. These results have larger uncertainties but there 
are no systematic differences from those shown in Figs 3 & 4. 

The main results of this work are summarized in Fig. 5 which shows how the solutions 
vary as a function £7a or degree of sub-Hubble expansion. The statistic in the individual 
panels in Fig. 5 was determined in exactly the same manner as that in the simulation 
in Fig. 2. Red points represent Y(H V ) and black points Y(H(z g ). To the left of the 
break at R/R v ~ 2 these ordinates are plotted against R(H V )/R V while to the right against 
R(H(z g )/ R v . The same groups were used to construct each panel. These included those 
with 2.6 < R v (h^Mpc) < 3.3 in the middle panel. The large scatter, particularly among 
the black points, is generally to values greater than unity. The most deviant among these 
points also tend to show the largest uncertainties. As discussed earlier, the full effects of 
the peculiar motions are assumed to be causing this bias. The red points represent a region 
with a cooler flow so that this bias should be negligible here. Giving most weight to the 
red points, we adopt the value of Vt\ ~ 0.7 and based on the spread of the mean of the red 
points from each panel assign an uncertainty of ±0.1. Using the more direct approach T08 
find Q\ ~ 0.77 from known distances and velocities. Finally we note from Fig. 1 of T10 
that the actual value of H v should increase slightly over the range 1 < R(H V )/R V < 1.7. 
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4. Summary and Conclusions 

From a combination of geometry and statistics we have shown that the velocity field 
surrounding groups of galaxies obeys two Hubble like expansion laws i.e. V = H exp R. At 
large distances the global Hubble Law is recovered. Within two well defined radii R v and 
Res we find a lower expansion rate (H exp ~ 0.8AH o ). Thus the main result of this work is 
the demonstration that the phenomenon reported by T08 extends well beyond the vicinity 
of the local group. If one adopts the the interpretation of T08 that the result is due to the 
local effects of dark energy, we find fl\ ~ (0.84) 2 ~ 0.7 ± 0.1 from a relatively small range 
of group mass. It is interesting that Sandage (1986), using the data then available, found 
that the calculated deceleration due to the Local Group is equivalent to H exp ~ 0.83H o at a 
distance ~ l.Qh^Mpc for an f2 = (A = 0) cosmology. 

These results do not contradict the conclusions from the cosmological simulations of 
Hoffman et al. (2008) & Martinez- Vaquero et al. (2009) that the 'coolness' of the local flow, 
by itself is not a strong discriminant between a model without A (OCDM) and ACDM. Here 
the flow between R v and Res is cooler because the galaxies within this region originated 
from inside the group and hence share the group's peculiar motion. Galaxies beyond Res feel 
the full effect of the peculiar motions. The observations show a systematic lowering of the 
expansion rate within the inner (cool flow) region in both of the OCDM and ACDM models. 
Based on the T08 interpretation the sub-Hubble rate in the ACDM model is predictable and 
when applied can account for the observations. However, interpreting the results as a local 
cosmological effect does contradict the conclusions from the above simulations. Until this 
impasse is definitively resolved with further work one may consider H v to be a convenient 
parameterization of a sub-Hubble expansion rate. 

In conclusion, significant perturbations to the global flow around galaxy groups beyond 
the Local Group are present. The sub-Hubble expansion rates are similar to those found 
recently by T08 around 3 nearby groups. The possibility that local dark energy effects 
are responsible makes the problem worthy of further investigation both observationally and 
theoretically. 
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Fig. 1. — The geometry of the problem is shown. The group is at distance d g while the field 

cralaw at r\ r tnihtpnrlt! an ancrlp cv at thp nVisprvpr Spp t.pvt. for Hptailt: 
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Fig. 2. — The results of a simulation showing the average V/H exp R versus normalized dis- 
tance from the group center R/R v . In this simulation group properties were the same as used 
in the analysis of the observations of the field galaxies. Here the 'field galaxies' were assigned 
random angles of 9, assigned random group-centric distances R and given appropriate line 
of sight redshifts due either to global Hubble or sub-Hubble expansion about these centers if 
R/R v was greater than or less than 2. Three hundred field galaxies were randomly assigned 
to each of the 12620 groups. The result shown is for an (Q m = 0.3, = 0.7, Qk — 0) 
cosmology. The red points then represent Y(H exp ) with H exp = H v = 0.S4H o while the black 
points show results with H exp = H(z g ). Error bars here and in all of the following diagrams 
represent the standard error of the mean. 
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Fig. 3. — Average V/H exp R plotted against R/R v for the cosmology and the assumed ex- 
pansion velocity indicated. Only groups with 2.6 < R v {h^Mpc) < 3.3 were included as it 
is within this interval that the break at R/R v ~ 2 was most distinct. This point represents 



-14- 



1.4 



K -i 2 

x 

a; 



1.2 < R/R < 2 



> 



1 

0.8 
0.6 



, Jim? xziiHhii&hi 




2 3 

R h^ 1 Mpc 



1.4 

x 
aj 

\ 1 

> 

0.8 
0.6 



2 < R/R < 5 



hh hl tjj\kl\l l 



3*" 



2 3 

R v h" 1 Mpc 



1.4 



K -i 2 

QJ 



> 



1 

0.8 
0.6 



1.2 < R/R v < 2 



is**' 

if 4-T i T i i I i i i i I i i i 



2 3 4 

R v h-J Mpc 



Fig. 4. — Top and middle panels: Average V/H exp R plotted against R v for all 12577 groups 
with R v < 5.7 and for the two regions of R/R v indicated. As in the other figures red 
points represent Y(H V ) and black points Y(H(z q ))- The bottom panel is similar to the top 
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Fig. 5. — The same as Fig.2 except now showing the results from an analysis of the group 
and field galaxy observations for three trial values of As in the other figures red points 
represent Y(H V ) and black points Y(H(z„)). Note that R v itself is dependent on the assumed 



